Propagation and control of nano-scale, magnetic droplet solitons 
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The propagation and controlled manipulation of strongly nonlinear, two-dimensional solitonic 
states in a thin, anisotropic ferromagnet is theoretically demonstrated. It has been recently pro- 
posed that spin polarized currents in a nanocontact device could be used to nucleate a stationary 
dissipative droplet soliton. Here, an external magnetic field is introduced to accelerate and control 
the propagation of the soliton in a damped medium. Soliton perturbation theory corroborated by 
two-dimensional micromagnetic simulations predicts several intriguing physical effects including the 
acceleration of a stationary soliton by a magnetic field gradient, the stabilization of a stationary 
droplet by a uniform control field in the absence of spin torque, and the ability to control the soli- 
ton's speed by use of a time varying, spatially uniform external field. Soliton propagation distances 
approach ten microns in low loss media suggesting that droplet solitons could be viable information 
carriers in future spintronic applications analogous to optical solitons in fiber optic communications. 



Nanomagnetism holds great promise for future spin 
based information storage and processing technologies 
PQ. One enabling physical effect is spin torque [5] which 
imparts angular momentum from a spin polarized cur- 
rent to a magnet. Spin torque forms the basis for tunable 
microwave nano-oscillators in confined nanopillar struc- 
tures [3] and nanocontacts abutting an extended ferro- 
magnet [1] . Most nanopillar dynamics can be reasonably 
described by single domain modeling .5] with the notable 
exception of gyrotropic vortex motion [6]. In contrast, 
nanocontacts enable the excitation of radiating [8] and 
localized [7] coherently precessing nonlinear wave states. 
The analysis of solitonic waves in nanocontact systems 
has predominantly been limited to either the weakly non- 
linear regime at threshold [9] or complex micromagnetic 
simulations [5J [TU] . 

It was recently proposed that a spin torque driven 
nanocontact could act as a soliton creator in a uni- 
axial ferromagnet with sufficiently strong perpendiular 
anistropy [11) . The resultant strongly nonlinear, coher- 
ently precessing state was termed a dissipative droplet 
soliton, the locally driven/uniformly damped cousin of 
the two-dimensional (2D), nontopological droplet soliton 
[12j . Prior numerical computations suggested that the 
conservative, stationary droplet could be generalized to 
a propagating solution [13) . This was verified by the con- 
struction of a stable, two-parameter family of propagat- 
ing droplet solutions for the conservative Landau-Lifshitz 
equation [2]. In this Letter, we use soliton perturba- 
tion theory to semi-analytically demonstrate the feasi- 
bility of sustaining, moving, and controlling a droplet 
soliton solely under the action of an external magnetic 
field. Modulation equations describing the evolution of 
the soliton's speed and precessional frequency in the pres- 
ence of damping and a temporally /spatially varying ex- 
ternal field are studied and the results are corroborated 
by 2D micromagnetic simulations. We show that a sta- 
tionary droplet can be accelerated by a field gradient. 



Once in motion, the soliton's speed can be controlled by 
a spatially uniform, time varying external field. A field 
gradient due to two nanowires can accelerate a soliton 
to propagate approximately ten microns in a low loss 
ferromagnet. Stationary droplets of any allowable fre- 
quency can be created from a sufficiently large, localized 
magnetic excitation, induced by a nanocontact or other- 
wise, and then stabilized by a linear feedback control field 
without the use of spin torque. This represents a new 
mechanism to study magnetic solitons without strong, 
spin torque induced perturbations. Analogous to opti- 
cal solitons in fiber optic telecommunications [15] . these 
results show that droplet solitons act as stable, control- 
lable, particle-like, precessing dipoles exhibiting intrigu- 
ing nonlinear physics and holding potential for spintronic 
applications. 

The model of magnetization dynamics we consider is 
the Landau-Lifshitz equation [5] 
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describing a thin, two-dimensional, unbounded ferromag- 
net with perpendicular, crystalline anisotropy sufficient 
to overcome the demagnetizing field so that the steady 
state is O = 0. The magnetization is characterized by its 
polar and azimuthal angles O and <I>, respectively, a > 
is the damping parameter, and ha — ho(x,t) is the mag- 
nitude of an external magnetic field pointing in the direc- 
tion normal to the film plane with positive values parallel 
to the magnetization steady state O = 0. Time, space, 
and fields are normalized by scaled versions of the Larmor 
frequency \"l\HoM s {Q — 1), exchange length L c ^/y/Q — 1, 
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and saturation magnetization M S (Q — 1), respectively, 
where Q = Hk/M s > 1 is the scaled anisotropy field 
or quality factor measuring the relative strength of the 
anisotropy energy to magnetostatic energy. For Co/Ni 
multilayer anisotropic ferromagnets used in recent exper- 
iments [TB]> the temporal scale and length scale are ap- 
proximately 27 ps and 17 nm, respectively, and a ss 0.01, 
Q R3 1.25, M s ss 650 kA/m. References to dimensional 
results use these parameter values. 

In what follows, we assume that a large amplitude, lo- 
calized excitation has been nucleated by a spin torque 
nanocontact |llj or some other means. The rest of this 
Letter is concerned with the manipulation of this struc- 
ture in a lossy medium using an external field. 

When ho = a = 0, eq. admits the conserva- 
tion of spin density TV — J (I — cosQ)dx, momentum 
V = J(cos6 - l)V$d£, and energy S = \ J[|V9| 2 + 
sin 2 @(1 + | V<I>| 2 )]<icc where all integrals are taken over the 
plane. Minimizing the energy subject to fixed TV and V 
leads to a two-parameter family of localized, precessing, 
stable traveling waves called propagating droplet solitons 
parameterized by their velocity V and frequency in the 
comoving frame uj 114] , There is a bijective map from 
(TV, V) to the physical parameters (w, V). Droplet local- 
ization requires that the velocity and frequency of the 
propagating droplet must lie below the spin wave band, 
enforcing the restriction 



+ \V | 2 /4 < 1. 
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Typical droplet widths are order one, hence are nanoscale 
excitations. Stationary droplets with rest frequencies 
close to zero resemble static circular bubbles which re- 
ceived a great deal of attention in the past [17] • However, 
typical bubble sizes are much larger. With the inclusion 
of nonlocal magnetostatic fields, it is predicted that a 
static bubble will be stable for a 5 nm thick film with 
a radius above 63 microns [15] . Thus, droplets can be 
viewed as smaller, dynamic generalizations of the static 
bubble. 

Allowing for weak damping (a <C 1) and a slowly 
varying magnetic field (] V/ioj, \dth \ <C 1) causes the 
spin density, momentum, and energy to evolve in time. 
Through the map to (w, V), we can describe the droplet's 
adiabatic, particle-like evolution by a time dependence of 
the droplet's rest frequency and velocity, trajectories in 
the V-uj phase plane. 

Using soliton perturbation theory (see, e.g., [15]), we 
obtain the following system of modulation equations de- 
scribing the slow modulation of the spin density and mo- 
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The energy £ = So + \ f ho(l — cos 9) dx is constrained 
to evolve according to 

% = (. + h o) d § + (O t ho + V-VhoW + V. d §, (4) 

thus it is sufficient to evolve eqs. ^ only. The integrals 
are evaluated with conservative droplets (9, $) of given 
spin density TV(i) and momentum V{t) or, equivalently, 
rest frequency uj(t) and velocity V(t). The field is evalu- 
ated along the soliton trajectory ho — ho(X(t),t) where 
dX/dt = V. Similar modulation equations were derived 
for one-dimensional droplets in [501 121] • Due to the ro- 
tational invariance of eq. ([T]), we limit further discussion 
to droplet motion in the x direction so that V = (V, 0) 
and V = (V, 0) without loss of generality. 

In |14j . we numerically compute a library of propa- 
gating droplets using an iterative technique [22]. These 
precomputed states are used to numerically solve the 
modulation equations ([3]) and to recover w(i), V(t). 
We also perform micromagnetic simulations of the vec- 
torial form of eqs. (JlJ using a pseudospectral method 
described in [T5]. To recover the micromagnetic solu- 
tion's speed and frequency, we compute the center of 
mass / x(l — cos 9) dx / 'TV and the phase at the center 
of mass at each time step. The velocity and comoving 
frequency are then found by differentiation. The rest 
frequency is recovered by subtracting the local magnetic 
field ho(X(t), t). The modulation eqs. ^ represent a low 
dimensional projection of the magnetodynamics enabling 
us to bring finite dimensional dynamical systems meth- 
ods and control theory to bear on the problem. As we 
demonstrate, the results from modulation theory agree 
exceptionally well with micromagnetics. 

The dynamics of eqs. ^ depend on the magnitude of 
\Vho\/a. When \Vho/a\ ^> 1, the spin density is approx- 
imately conserved and the momentum varies. The V-uj 
phase plane for ho = 0.5 — 0.005a; and a — 0.001 pictured 
in Fig.JTJa) closely resembles trajectories of constant spin 
density shown in [T3]. The stationary droplet is accel- 
erated while the rest frequency decreases, representing a 
transfer of the effective potential energy stored in the pre- 
cessional motion to to effective kinetic energy of transla- 
tional motion V. For the balance iV/io/al = 0(1), differ- 
ent dynamics occur. Figure [ljb) depicts trajectories for 
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FIG. 1. (a,b,c) Droplet trajectories in the V-ui plane from 
numerical integration of the modulation equations |3| (solid) 
and micromagnetics (dashed). The dotted curves correspond 
to the edges of the precomputed droplet library and the thick, 
solid parabolic curve is the spin wave band, (a) low loss 
a = 0.001. (b) a = 0.01, positive bias field, (c) a = 0.01, 
negative bias field. See text for further details, (d) Droplet 
velocity time dependence for the solid/dashed trajectories in 
(b,c). (e,f) Droplet profiles from micromagnetics correspond- 
ing to the circles in (b,c), respectively. The gray (color) scale 
represents the out of plane magnetization cosO and the ar- 
rows represent the in plane components. 



the same field as in [ija) but with a = 0.01. The dashed 
curve depicts the trajectory from micromagnetics with 
w(0) = 0.37. The circle on this trajectory corresponds 
to the droplet shown in Fig. [ije). The droplet is ac- 
celerated and is accompanied by an amplitude decrease 
until it devolves into a linear spin wave upon reaching 
the band edge (eq. j2j|). The micromagnetic simulation 
closely matches the adiabatic theory until the band edge 
is reached and the solution amplitude is very small, after 
which the droplet ansatz is no longer valid. A plot of 
V{t) is shown in Fig.[]Jd). 

During the course of evolution, a droplet can experi- 
ence deceleration as in Fig. [TJc) with a — 0.01 and a 
negative bias field, ho = —0.5 — 0.005a;. This behav- 



ior is reminiscent of Bloch oscillations predicted for one- 
dimensional droplets in 21J. In the ID case, the soli- 
ton oscillates under a constant force. We have not ob- 
served this behavior with our micromagnetic simulations 
nor modulation theory. Instead, micromagnetic simu- 
lations reveal the formation of local, topological struc- 
ture including vortex/anti-vortex pairs and what appears 
to be switching of the magnetization to the inverted, 
= 7r state, corresponding to (V, w) = (0, 0) as shown 
inFig.^f). 

In addition to the dynamics of Fig. [T] for constant field 
gradient, we have also simulated droplet acceleration due 
to the field generated by two current-carrying nanowires 
in the plane of the film with current in the same direction 
and a stationary droplet nucleated in between them |23j . 
For the moderate damping case, a = 0.01, the droplet is 
predicted to travel up to 3 /an in about 30 ns for 15 mA 
current in each wire with a 3.3 /im wire separation. Top 
speeds can approach 600 nm/ns. In the low loss case, 
a = 0.001, the droplet can propagate about 10 /im in 70 
ns for a 10 mA current with a 10.3 /im wire separation. 

The remaining regime when |V/io|/ a 1 is now in- 
vestigated. We focus on the case of a uniform and static 
magnetic field, assuming that a propagating droplet has 
been created. In Fig. [2] solution of the modulation sys- 
tem |3]) reveals the acceleration of a propagating droplet 
due to damping when the magnetic field is zero or posi- 
tive. When the field is sufficiently negative, the droplet 
can experience deceleration and then acceleration as its 
amplitude decays. This counter-intuitive droplet accel- 
eration due to damping was predicted for ID droplets 
in the absence of a magnetic field 20 . For negative 
fields, the magnet undergos a complete reversal to the 
(y,w) = (0,0) state for initial droplets with parameters 
lying below a separatrix (see Fig. [2jc)) . The switching 
separatrix corresponds to the stable manifold of the fixed 
point (V,uj) = (0, —ho). Linearization of eqs. ^ around 
this fixed point results in the eigenvalues 
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d v T 
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evaluated at the fixed point. Since d u N < [12] and 
dyV > for stationary droplets, the fixed point is a sad- 
dle. The switching separatrix from modulation theory 
accurately resolves the micromagnetic dynamics as evi- 
denced by the close agreement in Fig.[2jc) for trajectories 
starting very close to the separatrix. It is the existence 
of this switching separatrix for negative bias field that 
led to the differing phase plane trajectories in Figs, [ljb) 
and (c). By changing the external field, the switching 
separatrix moves, suggesting a mechanism to stabilize a 
stationary droplet at the saddle point (V,w) — (0,—ho) 
by use of feedback control. 

The nucleation of a stationary droplet by a spin torque 
driven nanocontact has been theoretically demonstrated 
[llj . The droplet, which would otherwise decay due 



FIG. 2. Droplet trajectories with uniform, static magnetic field and a — 0.01. (a) Modulation solution showing acceleration or 
deceleration of a propagating droplet due to damping and ho = —0.4 (solid), ho — (dashed), and ho = 0.9 (dash-dotted), (b) 
ho = 1. (c) ho = —0.7. In (b) and (c), the dashed curves are micromagnetic simulations. In (c), the dash-dotted curve is the 
separatrix between the switched state (V, cj) — (0, 0) and spin wave states. 



to damping, is sustained by a balance between local- 
ized driving and uniform damping. However, the cur- 
rent induced Oersted field strongly perturbs the station- 
ary droplet from its ideal, symmetric structure, lead- 
ing to phase variations and potentially a drift instabil- 
ity whereby the droplet is ejected from the nanocontact 
[11) . Furthermore, canting of the polarization layer is 
required to obtain an ac electrical signal via giant mag- 
netoresistance leading to symmetry breaking of the spin 
torque term and further complexity. We propose a sim- 
ple alternative stabilizing mechanism that avoids these 
difficulties: a closed-loop, spatially uniform control field. 
The stationary droplet saddle point (V, ui) = (0,w*) for 
ho = — is altered to an attractor by introducing the 
linear feedback control 

ho(t) = -u, + g[u* - w(t)], g > 1, (6) 

with gain g and bias — u;* . Stability analysis and numeri- 
cal computations demonstrate that (V, w) = (0, w«) with 
the control field Q is a global attractor [23] ■ We have 
also performed micromagnetic simulations incorporating 
the feedback control law §6§. Quite generic, localized, 
large amplitude initial conditions lead to relaxation to 
the stationary droplet with uj — Simulations showed 
droplet locking with a 2 ns control field response time 
[23] , This suggests that any sufficiently large, localized 
excitation created by spin torque or other means can be 
deformed into a stationary droplet of choice solely with 
the use of an external, spatially uniform magnetic field. 

In the presence of damping and a constant magnetic 
field, a propagating droplet is either accelerated or de- 
celerated (recall Fig. [2]) . A time varying, spatially uni- 
form open loop control field can be used to stabilize the 
droplet's speed. We compute the field required to main- 
tain constant speed in the modulation equations ^ and 
use this field profile in micromagnetic simulations. Re- 
markably, the control field from modulation theory leads 



to accurate control of the droplet's speed in the full mi- 
cromagnetic simulation [23]. 

We have demonstrated that droplets propagating in a 
realistic, damped ferromagnet can be sustained, acceler- 
ated, and controlled by use of only an external magnetic 
field. Their robustness and controllability hold promise 
for future spintronic applications. 
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ADDITIONAL MATERIAL FOR 
"PROPAGATION AND CONTROL OF 
NANO-SCALE, MAGNETIC DROPLET 
SOLITONS" 



ACCELERATING DROPLETS WITH 
NANOWIRES 



Figure [3] depicts the acceleration of a nucleated droplet 
by means of two in-plane nanowires. The current induced 
Oersted fields lead to a negative magnetic field gradient 
that accelerates the droplet. The modulation system, 
eq. 3, enables a much larger investigation of parameter 
space than micromagnetic simulations. The simulations 
incorporate the Oersted field due to two infinite wires 
with 150 nm diameters and varying separation. Station- 
ary droplets with varying frequencies are assumed to be 
nucleated 300 nm from the center of the left wire and al- 
lowed to propagate until either reaching the second wire 
or the soliton center of mass attains the value m z = 0.5. 



Oersted field 



current 




FIG. 3. Droplet acceleration by current flowing through two 



nanowires. 



STATIONARY DROPLET STABILIZATION 

A stationary droplet nucleated by a spin torque driven 
nanocontact experiences significant perturbations due to 
the Oersted field and canting of the applied field and 
fixed layer orientation that can lead to a drift instabil- 
ity [IT]. To avoid these difficulties, we propose a stabi- 
lization mechanism based on linear feedback that turns 
the stationary droplet saddle point (V,uj) = (0,0;*) for 
ho = — uj* into a stable fixed point via 



/*>(*) 



(6) 



This control could be implemented using a small am- 
plitude (sub-threshold) dc current applied to a trilayer 
nanocontact with a canted fixed layer resulting in a mea- 
surement of f2 = u + ho with negligible spin torque and 
Oersted field effects. The goal is to drive f2 to zero. With 
such a control law, linearization of eq. 3 around the fixed 
point results in the eigenvalues 
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Since d u Af < and d v V > [14 , when g > 1 both 
eigenvalues are negative and the fixed point is linearly 
stable. This proves linear (spectral) stability. 

We numerically observe that the fixed point is a global 
attractor in Fig. |4j Figure [4ja) shows the relaxation of 
the frequency to = 0.5 for stationary droplets. Figure 
(4^b) shows that phase perturbations leading to droplet 
propagation decay in time hence the drift instability has 
been removed. The gain g determines the relaxation rate 
to the fixed point with larger values leading to faster 
relaxation. The constraint h a > — 1 is required in order 
to avoid spontaneous magnetization reversal in the far 
field. This limits the choice of the gain to 



g(uj* — uj(t)) > -1 + cj„ 



(8) 



If V(0) = 0, and the control law ^ is assumed, then 
eq. (3) simplifies to 



duj (1 — g)(cj* — oj) 
~dt =a KjV 



sin 2 O dx, 
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differing bandwidths or update times so that ho(t) is up- 
dated periodically. Figure [5] depicts the evolution of the 

0.3 




FIG. 4. Stabilization of a stationary droplet with = 0.5 by 
the linear feedback control law |6| with a — 2. (a) Stationary 
droplet relaxation, (b) Trajectories in the V-u> plane. 



with V(t) = 0. Since d u Af < [12] and g > 1, we observe 
that w = uj* is a global attractor for fixed V(t) = with 
uj(t) relaxing monotonically to a;*. This implies that for 
w(0) < w*, inequality {[sj) will be satisfied for any g > 1. 
However, when w(0) > a;*, then one must choose g such 
that 1 < g < (1 — cj*)/(o>(0) — o>*). For small amplitude 
excitations with w(0) close to 1, g will be very close to 
1. Thus there is a trade-off between the relaxation time 
and the stability of the ferromagnet. 

We have also performed micromagnetic simulations of 
the vectorial form of the Landau-Lifshitz equation (1) in- 
corporating the feedback control law ([6]). We begin the 
computation with an asymmetric, localized initial condi- 
tion 

Q(x,y,0) = Aa-W^'-W"*? , $(x,y,0)=Q, (10) 

with A = 2.7, w x — 2.3, w y = 3. Equation (1) is evolved 
with the spatially uniform field Q where = 0.4, 
a = 0.01, 5 = 3. We "measure" = w(t) + /i (*) 

by averaging the in-plane magnetization orientation over 
the unit disk to obtain an average phase $(t). This mod- 
els the nanocontact measurement technique suggested. 
Then, fi(t) = ^§r{t)- We perform computations with 




2000 



FIG. 5. Micromagnetic simulation giving the spatially av- 
eraged magnetization frequency due to the control field |6}, 
uj, — 0.4, g = 3, a = 0.01, and an update time of 74. The 
droplet is locked when Q — 0. 



total frequency £l(t) with a field update time of 74 which 
translates to 2 ns or a 500 MHz bandwidth for perpen- 
dicular magnets used in recent experiments [16] . Because 
damping drives the dynamics, we expect that the opera- 
ble control bandwidth is approximately a\j\^ioM s (Q — 1), 
corresponding to about 400 MHz. This suggests that 
any sufficiently large, localized excitation created by spin 
torque or other means can be deformed into a stationary 
droplet of choice solely with the use of an external, spa- 
tially uniform magnetic field. 

Physically, this stabilization mechanism balances the 
switching of a droplet by a sufficiently negative bias field 
and the decay of a droplet via damping. 



DROPLET SPEED CONTROL 

In the presence of damping and a constant magnetic 
field, a propagating droplet is either accelerated or de- 
celerated (recall Fig. 2). A time varying, spatially uni- 
form open loop control field can be used to stabilize the 
droplet's speed. For this, we implement an optimization 
strategy for the modulation equations (3) by stepping for- 
ward in time and determining the field h (t) that enforces 
the constraint of constant speed. The resulting field pro- 
file is fit to a quintic polynomial (see Fig. |6^b)) and used 
in a micromagnetic simulation of eq. (1) resulting in the 
trajectories shown in Fig. |6^a). Remarkably, the control 
field from modulation theory leads to accurate control of 
the droplet's speed in the full micromagnetic simulation. 
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FIG. 6. Open loop control of droplet speed, (a) Two trajecto- 
ries (solid, dash-dotted) in the V-ui plane from micromagnetic 
simulations with (b) the corresponding control field (solid, 
dash-dotted, respectively) determined from the modulation 
equations (3). 




